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Abstract 

We study the free energy landscapes of a pair of submicron spherical particles floating at the 
surface of a sessile droplet. The particles are subjected to radial external forces resulting in a 
deformation of the droplet shape relative to the reference shape of a spherical cap. This deformation 
leads to tangential forces on the particles. For small deformations and for the contact angle 0q at 
the substrate being equal to vr/2, the corresponding linearized Young-Laplace equation is solved 
analytically. The solution is constructed by employing the method of images from electrostatics, 
where each of the particles plays the role of a capillary monopole and the substrate is replaced 
by a virtual drop with image charges and by imposing the conditions of fixed droplet volume and 
vanishing total force on the droplet. The substrate boundary conditions determine the signs of the 
image capillary charges and therefore also the strength of the tangential forces on the particles. In 
the cases of an arbitrary contact angle 9q these forces are calculated numerically by employing a 
finite element method to find the equilibrium shape of the droplet for those configurations in which 
the particles are close to the local free energy minima. 

PACS numbers: 68.03.Cd, 47.85.-g, 89.90. +n, 83.80.Hj 
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I. INTRODUCTION 



Present day chemical synthesis methods allow for the creation of colloidal particles of 
various sizes, types, and shapes The behavior of these tailored particles in suspensions 
mimics that of atoms at favorable length and time scales, thus providing model systems for 
testingbasic theoretical concepts, which even permit direct visualizations of thermal fluctua- 
tions [2|. In addition colloidal particles offer very rich perspectives for potential application, 
e.g., in the context of optically active materials [3j. Colloidal particles of intermediate wetta- 
bility usually adsorb at fluid-fluid interfaces. This configuration is extremely stable against 
thermal fluctuations due to the very large free energy of adsorption for micron-size particles. 
This can be used, for example, to stabilize emulsions J|. Such effectively two-dimensional 
colloidal systems have also been used to successfully test ^| the Kosterlitz-Thouless j^] sce- 
nario for 2D-melting, which is very difficult to capture in atomic systems. The presence 
of an interface affects effective interactions between colloidal particles considerably as com- 
pared to those for colloids in bulk solution. Examples for these interface induced effective 
interactions include capillary forces or effective dipole-dipole repulsive forces at water-air 
interfaces Whereas the gravity induced former ones do not play a role for colloids with 
a size below a micron, the latter ones explain the emergence of 2D-crystal structures 8J for 
laterally confined systems. However, there are experimental evidences js-ll] supporting 
the occurrence of attractive forces between equally charged particles at interfaces, with a 
range exceeding by far the range of van der Waals attraction. The observed attraction was 
originally attributed to flotation-like capillary forces decaying proportional to d~ x , where 



d is the distance_between the particles 
unsatisfactory 



111 ] . However, this explanation was found to be 



12], because for a mechanically isolated system, composed of the particles 



and the interface, one rather obtains an attractive behavior oc d~ A [13|. Consequently, 
the effective capillary attraction cannot overcome the direct dipole-dipole repulsion oc d~ A , 
unless the screening length in water is comparable with the size of the particles. In this case 
a total effective interaction potential exhibits a shallow minimum with a depth of several 



131 ] . In Ref. 12| it has been suggested that the finite curvature and the pinning of 



the interface might be important for generating a long-ranged attraction oc d 1 . In fact, in 



the experiment reported in Ref. 



Ill ] the particles were trapped at the surface of a droplet 



pinned to a solid plate, but this aspect was neglected in the theoretical analysis put forward 
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in Ref. [11]. 

Recently, it has been shown [l^] that finite curvature and pinning indeed have profound 
effects on the trapping potential of a colloid as compared to the case of a flat and unbounded 
interface. In addition, general expressions have been obtained for the free energy of the 
system in the case of an arbitrary external pressure field acting on the droplet surface. In 
the present analysis we use those results in order to calculate the effective pair potential for 
particles subjected to external radial forces (see Fig. 

We consider N spherical particles floating at the surface of a sessile, non- volatile droplet, 
the surface of which meets the planar substrate at a contact angle 6*0. With no external forces 
acting on the particles the shape of the droplet remains that of an undeformed spherical 
cap of radius Rq and the immersions of the particles into the liquid forming the droplet are 
determined by the contact angles at each of the particle surfaces 151 ] . (Due to the partial 
immersion of the colloidal particles and the volume constraint of the liquid, the radius Rq 
of the spherical cap with colloids is larger than that of the corresponding droplet without 
colloids.) We call this the reference configuration. We place the origin of the reference 
system at the geometrical center of a sphere formed by completing the unperturbed droplet 
shape to a full sphere, and choose the outward normal of the exposed substrate surface as the 
positive z— direction. We parametrize the positions of the particles by the radial distances T{ 
and by the solid angles fi, = (6^, (pi), where Q = (8, <fi) is the solid angle used to parameterize 
the droplet surface (see Fig. [T|). In the reference configuration the particles float freely at the 
surface of the droplet, because in this case there is no preferred angular configuration of the 
particles. In the presence of radial external forces fi& r (e.g., generated by optical tweezers) 
acting on each of the particles the droplet deforms such that those forces are balanced by the 
radial components of the corresponding capillary forces. However, the tangential capillary 
forces acting on the particles are in general unbalanced which leads to effective interactions 
between the particles. In terms of the free energy this means that in this case the total free 
energy depends on the angular configuration of the particles. 
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FIG. 1: (Color online) Geometry of the system under consideration. The positions of the colloidal 
particles trapped at the surface of a sessile droplet with contact angle 9q are parameterized by 
solid angles f^i and f^- These particles are subjected to radial external forces fi = /ie r (f2i) 
and f*2 = fz&ri^), respectively, where e r is the unit vector normal to the spherical cap-like 
surface of the droplet in the reference configuration. In Ref. [ljj it has been shown that under 
these conditions the particles can be well approximated by capillary monopoles, which in turn are 
determined exclusively by the external forces. This implies that the sizes of the particles and the 
contact angles at their surfaces are irrelevant. The dashed lines on the surface of the cap of the 
sphere with origin O and radius Ro are the points with spherical coordinates (Ro,0i,(j>i), i= 1,2, 
with (pi fixed. #21 is the angle formed by the vectors f*i and ¥2- The angle 021 is defined in that 
tangential plane to the spherical cap for which fi is a normal vector. At position 1 two curves 
meet: one is the intersection of the plane 4>\ = const through the origin O with the spherical cap 
(dashed line), the other is the intersection of the spherical cap with the plane through O spanned 
by fi and f*2 (long red arrow); ^21 is the angle formed by the vectors tangential to these two curves. 
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II. FORMULATION OF THE MODEL IN TERMS OF CAPILLARY 
MONOPOLES, GREEN'S FUNCTIONS, AND THE FREE ENERGY 

Previously it has been shown [l^] that the influence of a single colloidal particle on the 
surface of a sessile droplet can be described in terms of a capillary monopole determined by 
the external radial force acting on the particle. The results of this approximate, linearized 
theory agree very well with the corresponding results of the numerical minimization of the 
full, non-Gaussian surface free energy. Therefore in the following we describe the particles in 
terms of capillary monopoles, so that the particle sizes and contact angles at their surfaces 
remain intrinsic parameters which do not enter the analysis explicitly. Accordingly, we study 
the droplet under the action of radial external pointlike forces of magnitudes fi giving rise 

to an external pressure U(Q) = J^ti /i<*( n > n i)/ R o with W = 6 ( e ~ Q')K<t> ~ </>')/ 
as Dirac's angular delta function. 

We consider the case that the deformations u(Q) = R(Q) — Rq of the droplet surface R(Q) 
due to the external pressure IT(f2) are small. For this to hold, the condition [l^ <C 7-Ro, 
i — 1, . . . , N, where 7 is the tension of the droplet surface, might not be sufficient, because 
the forces on individual particles, even if they are small, might add up to a large total 
external force (and thus might lead to a large deformation); therefore we impose the more 
stringent condition Yli=i \fi\ ^ 7-Ro- If this condition holds, e := ^ i=1 \fi\/(jRo) can be 
regarded as a small parameter. (In the one-particle case this reduces to e = //(7-R0) as 



introduced in Ref. 



14j|). On the other hand one has to keep in mind that the forces /j 



cannot be stronger than the maximal capillary forces, which lead to the extraction of the 
individual particles from the surface, and which are of the order of 70$, where is the 
characteristic size of the particle i. 

The resulting shape of the droplet can be determined by minimizing a free energy func- 
tional under the constraint of a fixed liquid volume, which leads to the Young-Laplace 
equation relating the curvature of the droplet surface to the pressure difference across the 
interface [yj]. This equation depends also on the mechanism fixing the lateral position of 
the droplet at the substrate. In the following we shall use the subscript a in order to dis- 
tinguish between the cases of a free contact line with a fixed center of mass (a = A) and 
of a pinned contact line with a free center of mass (er = B; in this case the center of mass 
can be free because the balance of forces is already guaranteed by the pinning of the contact 
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line). Linearized in e the Young-Laplace equation, expressed in terms of the dimensionless 
deformation v := lim^o u/(R e), takes the form (see Eq. (29) in Ref. jl^t ) 

- (V* + 2)v(n) = 7r(fi) + 7r CM (fi) + fi, (1) 

g- 

where V a := egde H -<9</, is the dimensionless angular gradient on the unit sphere and u 

sin 

is the spatially uniform dimensionless shift of the internal droplet pressure relative to the 
Laplace pressure 2'j/R of a spherical droplet. This shift occurs due to the external excess 
surface pressure 7r(f2) := lim e _^ U(Q)R / (^ye) exerted by the external forces /j. The free 
energy approach reveals that the internal pressure is equal to A, where —A is the Lagrange 
multiplier introduced into the free energy functional in order to fix the liquid volume. There- 
fore \x can be identified as fi = lim e _ i .o(A(e) — 2 r y/R )R /('je). The effective pressure ttcm(Q) 
corresponds to a body force fixing the center of mass of the droplet, which has to be intro- 
duced in the case of a free contact line (a = A) in order to achieve mechanical equilibrium; 
accordingly, it depends on the total external force YliLi fi e i(&-i)- m the case °f a f ree center 
of mass (cr = B) one has ttcm — 0. Equation ([1]) is supplemented by the boundary condition 
sin OqOqv \g Q — cos6 v 1 O = for the case of a free contact line and v\g = for the case of a 
pinned contact line (Eqs. (32) and (33) in Ref. [l^]) and by the incompressibility condition 

dnv(n) = 0, (2) 

where the angular domain Q corresponds to the shape of the reference droplet, which is a 
spherical cap with the contact angle 9q at the substrate. 
The formal solution v = v a of Eq. (CQ) can be written as 

v a (n)= [ dn'u(n')G a (n,n',e ), (3) 

Ja 

where Green's functions G a satisfy 

- (v* + 2)G a (Q, n', o ) = 5(Q, n') + A a (n, n', e ). (4) 

The functions A a (Q, Q', 6q) correspond to /i and hqm in Eq. (CQ), so that /i + ticm(Q) = 
J Qo dQ' 7r(n')A .(n, Q' , 6$)_. They can be determined from the force balance and the incom- 



pressibility condition 14J. Here, we only point out that, in general, these functions are 



not symmetric with respect to interchanging Q and fl ; and so neither are Green's functions 
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G a (Q, Q', 9q). The dependence on a enters via the boundary conditions for G a : 

sin 9 deG A (tt, 9 )\aedn - cos 9 G A (Q, 9 ) \nedn a = 0, (5) 

G B (fi,n',0o)|nean o = O. (6) 

The incompressibility condition can be expressed as 

dnG a (n,n',9 Q ) = 0, (7) 



'n 

see Eqs. (50)-(52) in Ref. Q. We point out that, locally, i.e., for -)• fi', Eq. (JH) can be 
approximated as — V„G CT (f2, 6*o) = 5(f2, from which it follows that G a (Q, Q', 9q) > 



-ln(#)/(2vr) where 9 is t 
tt' (see Eq. (57) in Ref. 



re angle between the unit vectors pointing into the directions Q and 



16 



17j, this 



By adopting the analogy with 2D electrostatics 
logarithmic divergence can be traced back to the diverging self-energy F se if of a monopole 
associated with a pointlike force. However, in an actual system the size a of the particle 
serves as a natural cutoff for the divergence of the self-energy which for a — > scales as 
F self ~ f 2 \n(R /a) < ^ 2 a 2 hx{R Q /a), using / < 7 a. 

The excess free energy associated with the droplet deformation due to the external pres- 
sure 7r(f2) can be written in the form (compare with Eq. (53) in Ref. [jjj]) 

F a = --eV?fj f dtt [ dn'Tr{n)G a {n,n',9 )ir{Q'). (8) 
Inserting n(Q) = J2iLi fi^fo, ^WEili \fi\) one obtains 

N N 

i=l i=l i<j 

where F^ se if = f 2 /(4ii r y) ln(i? /aj) + 0(1) is the self-energy of particle i, which depends 
neither on the position of the particle on the droplet nor, in leading order in ai/R , on the 
contact angle 9q and the boundary conditions a. The one-particle free energy landscapes 
Ai 7 ^ are given by 

AF^(9 U 9 ) = -^[G a . reg (n i ,n u 9 )-G a , reg (0,0,9o)} = -|V(Mo) - 9a(0,9 )] : (10) 

where the function g a is independent of i and corresponds to the regular part G atTeg (i.e., not 
containing the logarithmic divergence) of Green's function which can be found analytically 
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in the case 9q = tt/2 (see Sec. IV. C in Ref. [l^j and Fig. [2]) and numerically otherwise 
The last term in Eq. consists of effective pair potentials V a given by 

VriSh, Qj, 9 ) = -^0 [G CT (a, Qj, 9 ) + G^fy, n t , 9 )] . (11) 



We note that V a is explicitly symmetric with respect to fli and flj even if G a is not. 

For the special case 9n = tt/2 one can construct G a by using the method of images from 

n n 

electrostatics [14j (see Sec. IV in Ref. [14J). In the reference configuration the substrate 
is replaced by a virtual mirror image of the actual droplet, such that the union of the 
actual drop and of the virtual drop forms a full spherical droplet. The deformations of a 



17], who derived the following 



fully spherical drop have been studied by Morse and Witten 
equation for a "free" Green's function G: 

- iyl + 2)G(n, no = £ Y^WJp) = 5(n, n'), (12) 

l>2,m 

where the right hand side is an expansion in terms of spherical harmonics of a modified 
Dirac delta function 5(Q, Q') with the I = and I = 1 components projected out. Summing 
up the coefficients in the expansion in terms of spherical harmonics one obtains the solution 
in closed form: 



G (n,n') = G(9) = -^ 



14- - A-cos0 
— I — cos v + cos v In 



(13) 



2 3 V 2 

In the case of a sessile drop the solution is constructed in the form of a superposition of 
the point-force solution given by Eq. ( II 3p (which would be valid, if there was no substrate) 
and a solution associated with a virtual pointlike force acting symmetrically at the surface 
of the virtual lower hemisphere. For constructing Green's functions one has to take into 
account the boundary conditions as well as the mechanisms imposing the corresponding 
force balance. In the case of a free contact line the droplet is not attached to the substrate 
and mechanical equilibrium is achieved by fixing the lateral position of the center of mass of 
the liquid. This latter aspect is already incorporated into Green's function as given by Eq. 
(I12p . However, if the balance of forces is imposed by a pinned contact line and the center 
of mass is not fixed, the corresponding theoretical treatment requires the introduction of 
additional images, so that the total force vanishes. Concerning these details we refer the 



reader to Refs. 



14) and [17]. 
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III. RESULTS 



A. Contact angle 6q = 7r /2 
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FIG. 2: The functions —g a (0i, 9q = vr/2) +g a (0, 6q = vr/2), a = A,B, which determine the behavior 
of the one-particle free energy landscape AF^ (6i,9o) of particle i, as functions of its polar angle 
Oi, for 6 = vr/2 (see Eq. (fTU|)). 



First, we present our analytical results for the case 8 = n/2. As can be inferred from 
Eq. ( ITOj) and Fig. |2] the one-particle free energy landscape AF^ (8i,8o = 7r/2) of particle 
i is a non-monotonic function of its polar angle 6*j. As a consequence, besides the known 
phenomena of attraction of a particle to a free contact line (model A) and repulsion from a 
pinned one (model B), one finds a local free energy minimum for the particle being at the 
drop apex Oi = and a local maximum at 8{ = 6 max w 52° for model A, and a minimum 
at 8i = Q m in ~ 49° for model S. In the case of many particles the one-particle free energy 
landscapes compete with the effective pair potentials V a . The quantity relevant for obtaining 
the actual configurations of the particles is the excess free energy defined as 



N 



N 



AF W :=F W_£ Fi 

i=l 



/ - i,self 



(14) 



i=l 



i<j 



where we have introduced the notation V a ^ = V a (Qi,flj,9 ). In Figs. [3] and |4] we present 
the results for AF& in the case N = 2. We consider particles placed at Q\ = 
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and f2 2 = (# 2 , 02) and subjected to external radial forces f\ and / 2 such that |/i| <C 7-Ro, 
I/2 1 <S 7-Ro, and |/i + /2I <C 7-Ro- The excess free energy AFg is calculated for a fixed 
angular position of the first particle, referred to as the reference particle, as a function 
of the angular position of the second one acting as a probe particle. The angular cutoff 
S determines the closest approach of the particles to each other, i.e., #21 > 5, and to the 
contact line, i.e., 9i, 62 € [0, 0q — 5], where #21 is the angle between the vectors pointing in the 
directions Q\ and Q2 (see Fig. [TJ). One has 5 > 2a/ Rq, but one has to remember that for the 
configurations corresponding to #21 — 5 one should expect that the actual behavior deviates 
from that obtained within the monopole approximation. In the following calculations we 
take 5 = 7r/36, which corresponds to Ro ~ 23a. 

The interaction potential V a is a function of the angular coordinates of both particles 
separately, i.e., it is not only a function of their separation. This is important for all angular 
configurations due to the long range of the capillary deformation around a monopole. This 
deformation does not depend on the radius Rq of the droplet but only on the strength 
of the external force. In those cases in which the particles would correspond to higher 
capillary multipoles the pair interactions would vanish for Rq — > 00; in this sense the case of 
monopoles is exceptional. We also note that, as given by Eq. ffTTj) . V a is explicitly symmetric 
with respect to Vt\ and Q2, such that the aforementioned asymmetry of Green's functions 
G a is not proliferated to the free energy. 

Due to the rotational symmetry of the reference droplet forming a spherical cap the free 

(2) (2) 

energy depends only on the difference 02 ~ 0i so that AFo- = AFr(0i,0 2 ,0i - 2 ). We 
also introduce an auxiliary azimuthal angle 02i which describes the angular position of the 
probe particle relative to the reference one (see Fig. [1]) and thus the orientation of the pair 
of particles on the droplet. 

In the case of a free contact line (Fig. [3]) at the substrate and Q\ = the minima of the free 
energy are degenerate both for the configurations with the probe particle at 62 = 5 (i.e., in 
contact with the reference particle) and for those with 82 = 7r/2 — 5 (i.e., the probe particle 
being in contact with the substrate). In both cases the dependence on 02 is degenerate. 
However, if Q\ ^ 0, the rotational symmetry is broken and this degeneracy is lifted. The 
precise values of the angular positions of the free energy minima depend on the value of 81 
but in general they correspond either to the closest approach of the particles to each other or 
to the contact line. If Q\ <§C 1, there are only two free energy minima. They occur at 6*21 = 5, 
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FIG. 3: Color coded effective pair potential for the case 9q = vr/2 and for a free contact line at the 
substrate (a = A). The colors represent the excess surface free energy 7AF4 /(/1/2) (Eqs. ([9]) 
and (|14p ) of the system as a function of the position VI2 of the probe particle for a fixed position 
ill °f the reference particle (white cross). Each pair of panels in a column corresponds to the 
same configuration. The white dashed lines schematically indicate those positions Q 2 of the probe 
particle which are at a fixed minimal polar angular separation 621 = 5 = 7r/36 from the reference 
particle, i.e., at mutual contact. Along these dashed lines Ai 7 ^ is locally minimal (maximal) at 
the black (white) dots. The black symbols (#,I,A) correspond to minima under the constraints 
621 > 5 and 2 < n/2 - 5. (a), fej[back side]: = (6>i = tt/18>i = 0); (b), f/J[back side]: 
n t = (6»i = 3tt/18,0i = 0); (c), ^[back side]: = (6>i = 5tt/18,^i = 0); (d), ^[back side]: 
^1 = (#1 = 77r/18, 4>i =0). Note that in (d) two minima are shown (#,A) separated by a potential 
ridge. In (c) the value 0\ = 57r/18 = tt/3.6 is taken slightly smaller than Q m ax ~ vr/3.46, which 
is the position of the free energy maximum for the reference particle alone (see Fig. [2]). In such 
a case the upper black dot marks the deeper one of the two local minima on the dashed curve. 
The two local maxima (white dots) are equally high. The orientations of the coordinate axes in 
(b)-(d)[(f)-(h)] are the same as in (a)[(e)]. The quoted numbers of minima are the total ones 
occurring for the given configurations shown in the columns of panels. 
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FIG. 4: Same as Fig. [3] for a pinned contact line at the substrate (a = B). (a), (d)\b&ck side]: 
fii = (6»i = 3tt/18,0i = 0); f&j, fej[back side]: Oi = (6>i = 5vr/18,0i = 0); fcj, f/J[back side]: 
^i = (#i = 7tt/18, 0i = 0). In ^ the value #i = 57r/18 = 7r/3.6 is taken almost equal to 
&min ~ vr/3.67 which is the position of the free energy minimum for the reference particle alone 
(see Fig. [2]). The positions J7i (white crosses) in (a), (b), and (c) equal those in Figs, (c), 
and (d), respectively. In (b) the two minima on the dashed white curve are degenerate. The lower 
white dot marks the higher local maximum on the dashed white curve. 

02i = 7r [•] (see 6*i = 7r/6 in Fig. Eta)) and at 6> 2 = n/2 — 5, 2 — 0i = ^ ffl (see 6\ = n/18 
in Fig. [6]^a)). Increasing 8\ leads first to the emergence of a third minimum at # 2 = vr/2 — 5, 
02 — 0i = [A] (see Q\ = n/6 in Fig. [6](a)) followed by the emergence of a fourth one at 
9 2 i = 5, 02i = [•] (see 6»i = 9 max « 52° in Fig. [5](a); for 9 max compare Fig. [2]). Finally, 
upon further increasing 8 1 the minimum at 9 2 i = 5, 2 i = n turns into a local maximum, 
leaving three local minima at (#21 = S, 2 i =0) [•], (# 2 = 7r/2 — 5, 2 — 0i = 0) [A], and at 
(62 = tt/2 — 5, 02 — 0i = 7r) [■] (see #1 = tt/3 in Figs. [5]^a) and[6](a)). The answer to the 
question which configuration corresponds to the global free energy minimum with respect to 
the positions of both particles depends on 5. For 5 <C 1 it is the configuration in which both 
particles touch the substrate and touch each other (see 61 = 47r/9, 2 — 0i = in Fig. [6]^a); 
for 6q = 7r/2 one has 6i < 7r/2), because the depth of the corresponding minimum increases 
as ~ ln(l/o~). Thus the particles arrange themselves parallel to the contact line. 

For a pinned contact line (Fig. H]) at the substrate we also observe a degenerate minimum 
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(a) free contact line (a = A), 9 21 = 5 = 71/36 



-0.25 f 



.cm -0.3 



: = 71/6 
1 = 71/4 

e , » e - 7C/3.46 

1 max 



e x = 571/18 

1 = 7C/3 
0^771/18 








3 4 
*21 



(b) pinned contact line (o = B), 21 = 8 = tc/36 
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FIG. 5: Normalized free energy AF^ (Eqs. (HD-fHH) and ([14])) of a pair of colloidal particles 
floating on a sessile droplet with contact angle 9q = 7r /2 and exposed to external radial forces fi 
as a function of the angle 021 (see Fig. [1]) for several fixed values of 9\ characterizing the polar 
angular position of the reference particle. Concerning 9 max and Q m in see Fig. [2j The polar angular 
separation #21 (see Fig. [1]) between two particles is kept constant upon varying <j>2\ and equals 
#21 = 6 = 7r/36, i.e., the free energy landscape is probed around the reference particle fixed in 
space, (a) free contact line; (b) pinned contact line. 
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(a) free contact line (o = A) 
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FIG. 6: Normalized free energy (Eqs. ([8l)- (llip and {H}) of a pair of colloidal particles floating 
on a sessile droplet with contact angle 9q = 7r /2 and exposed to external radial forces /, as a 
function of the relative azimuthal angle fa — <j>i (see Fig. [T])) for several fixed polar angles 9\ 
of the reference particle. Here, the free energy landscape is probed close to the substrate, i.e., 
O2 = Oq — 5 = tt/2 — ir/ 36. (a) free contact line; (6) pinned contact line. 
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of the free energy for 6± = and a broken symmetry for Q\ ^ 0. Upon further increasing 9\, 
first there is only a single free energy minimum, occurring at #21 = $, 4>2i = (see Q\ — tt/6 
in Fig. EJfc)). Upon increasing 9\ further this minimum splits continuously into two minima 
at (6»2i = S,(j)' 21 ) with <p' 21 e [0,7r] and at (6 21 = S,(j) 21 ) with <p 21 G [7r,27r] (see Fig. \§[b) for 
1 > 7r/4). For 6 1 ! = # min the minima have reached the values <p' 21 ~ n/2 and ^ ~ 37r/2, 
which corresponds to the configuration in which the particles are positioned parallel to the 
contact line (see Q\ = Q m in ~ 49° in Fig. 0(6)). Finally, for even larger Q\ the two minima 
merge into a single minimum at #21 = 5, 4>2i = ?r (see 9\ = n/3 in Fig. G3(6)). The capillary 
forces repel the particles from the contact line (see Fig. [2]) and therefore all configurations 
with any of the particles close to the contact line are energetically unfavorable. However, 
it might happen that one of the particles gets trapped close to the contact line by other 
means, for example due to an evaporative flux of the liquid towards the contact line or by 
adhesion to the substrate. If the probe particle be the trapped one, such that 9 2 = n/2 — 5 
and 02 is free, and for any fixed polar position 9\ of the reference particle (see Fig. E(6)), the 
preferred position of the probe particle at the contact line always corresponds to <p2 — <P\ = 0, 
i.e., the two particles are positioned at a great circle perpendicular to the contact line. This 
means that for any angular position of the white cross in Fig. H] the minimum of the free 
energy along the contact line occurs at the point closest to the white cross. In the case that 
both particles are constrained to lie in the neighborhood of the contact line one observes a 
monotonic attraction (Fig. E(6), 9\ = 4tt/9), contrary to the case of a free contact line (Fig. 
HHa), 9\ = Att/9). The configuration corresponding to the global free energy minimum with 
respect to the positions of both particles without any constraints is such that the particles 
touch each other at d\ = 9 2 m 6 min pa 49° (see d\ = 9 min , 2 i = 7r ± 7r/2 in Fig. Ef6)). Thus 
the particles spontaneously arrange themselves parallel to the contact line at the common 
characteristic polar angle 9 min . 

In summary, in both cases the particles attract each other and, as a doublet, arrange 
themselves such that they are both placed as close as possible to the minimum of the one- 
particle trapping potential AF^ which, in the case of a free contact line, occurs at the 
apex and at the contact line whereas in the case of a pinned contact line it occurs at an 
intermediate angle Q m in- Additionally, in the case of a free contact line there is another local 
free energy minimum corresponding to both particles being at the contact line, however not 
touching each other but being positioned on the opposite sides of the droplet (<f> 2 — 4>i = 7r). 
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B. Arbitrary contact angle 



In the case of an arbitrary contact angle 8q 7^ tt /2 we minimize the free energy numerically 
by using a finite element method 



First, the one-particle free energy landscapes AF^-(8i, 8 ) are calculated for various values 
of the contact angle #0 and the force /,; acting on particle i. The numerical procedure is 



similar to the one described in Ref. 



14J and consists of a pre-evolution of a body of liquid 



towards the reference configuration (/j = 0) for a fixed polar angle 8i of particle i, after 
which the force is turned on (/j 7^ 0) and the droplet surface evolves towards the minimum 
energy while the particle is allowed to move only radially. We recover the scaling of the 
free energy ~ ff, which yields the scaling functions g a {.8i,8 ) (Eq. ffTUj) ). In the case of a 
free contact line we observe small deviations from this scaling, which we expect to vanish 
for very large droplets (which, however, are beyond the reach of the finite element method). 
There are two sources of this deviation. First, there is a finite-size effect due to fixing the 
center of mass of the droplet which leads to a contribution to the free energy associated with 
the work done by the virtual counterbalancing force fcM,i{@i) — /isin#j in displacing the 
center of mass from its reference lateral position x = xcM,ref(8i), which does not depend 
on fi, to x = (for details see Ref. Ijj]). Thus, this correction is linear in fa and it can be 
eliminated by taking the average of the results for f] = and fi = —\fi\- Second, there is a 
numerical error associated with the possibility that the evolving surfaces can get trapped in 
local free energy minima and thus do not reach the global free energy minimum. Performing 
the calculations for positive and negative values of fi and subsequent averaging minimizes 
this error, too. 

Qualitatively, the results (see Fig. [7]) are the same as in the case #0 = 7r/2 (Fig. [2]), i.e., for 
a free contact line the free energy Afj 1 ]^,,^) has two minima: at the drop apex (9i = 0) 
and at the contact line (6*,; = 8 — 5, where 5 is the closest possible angular approach of the 
particle to the contact line), whereas in the case of a pinned contact line AF^\(8i,8o) has 
only a single minimum at an intermediate angle 9{ = 6 m i n (8o), which strongly depends on 8q. 
Quantitatively, the energy barriers associated with the local free energy minima in both cases 
grow with 8q approaching n. According to our analytical theory, the radial displacement 
of the particle h is proportional to Green's function (see Eqs. (57) and (58) in Ref. [l^|): 
h = |/Ma)/7 + 0(a) = f[HRo/a) + 2<icG reg , ( ,(n i ,Sl i ,do) + 0(1)]/(2tt 7 ) (see Eq. ©). 



1(3 



Therefore the effective elastic modulus of the droplet, defined as k e ff := df / dh = (dh/df) 1 , 
is given by 

keff = k eff , a = 27r 7 /[ln(i?o/a) + 2tt^(^, 9 ) + 0(1)] 

= 27T 7 /{ln( J Ro/a) + 2ng a (0, 9 ) - 2n[-g a (9 i , 9 ) + g a (0, 9 )] + 0(1)]}. (15) 

Due to the good agreement of the analytical and numerical results for 9 = it/2 and for 
droplet sizes down to Ro/a — 4, as discussed in Ref. [ijj], one can expect that the corrections 
0(1) indicated in Eq. ffT5]) . due to the finite size of the particle, are small. Accordingly, our 
numerical results in Fig. [7] lead to the conclusion that in the case of a free contact line at 
the substrate (a = A) the droplet is relatively softest if the pulling force is either directed 
along the z-axis or applied at the contact line, whereas in the case of a pinned contact line 
(cr = B) the droplet is softest if the force is directed along the specific intermediate polar 
angle 9 min (9 ). 

In analogy to the case 9q = 7r /2, for two particles we expect that in order to reach 
equilibrium each particle drifts towards the polar angle 9{ close to that of the minimum of 
AF^(9i, 9 ). In the case of a free contact line this corresponds to a configuration in which 
both particles are positioned close to the drop apex or at the contact line. In the latter case 
we have evaluated the free energy numerically as a function of the azimuthal angle <p 2 — <fti 
(see Fig. [8](a)) for 9 X = 9 2 = 9 — 5. We have performed the calculations for #0 < n/ 2 with 
fi — I2 > and for 9q > n/2 with fx — f 2 < 0, where the restriction on the sign of the force 
follows from the constraint that the particles cannot penetrate the substrate. Having for a 
given #0 the data accessible only for one sign of the force we could not take the average; 
instead we had to apply a finite size correction (see Appendix A), analogous to the one 
proposed in Ref. [jjl]. The data show that there are one or two minima depending on 9$. 
There is always a minimum for <f) 2 — <fii = 5. If the contact angle #0 is smaller than a critical 
angle ps 5ir/6, there is a second minimum for 2 — <Pi = 71 '■ The free energy barrier between 
the two minima is largest for #0 ~ 7r/3 and it vanishes for 9q > 5tt/6. In the case of a pinned 
contact line (see Fig. M^b)) the free energy for the particles positioned at the polar angles 
9\ = 9 2 = 9 min (9Q) is, independently of #0, always a monotonic function of <p 2 — <j)\ with a 
minimum corresponding to the particles touching each other (<j) 2 — <pi = 5). The absence of a 
second minimum at <p 2 — <f>\ = ir in this case is rather surprising. For example, for #0 = 5tt/Q 
one has 9i = 9 2 = 9 min (9 = 5n/6) ~ ir/2, so that for 4> 2 — (j)i = vr the angular separation 9 21 

17 



(a) free contact line (a = 



O°0.3 
o 



+ 



0.2 



- e = 7i/6 

- e = 27i/6 

- 9 Q = 37t/6, theory 

- 9 = 47T/6 

- = 57c/6 



CD 



0.1 



CD 




(b) pinned contact line (a = B) 



CD 
o 

+ 



e = 7i/6 

9 = 2tc/6 



— 9 Q = 37C/6, theory 




FIG. 7: The numerical results for the scaling functions gcr(0i, 6q) (Eq. (I10p ) obtained by taking the 
average of the results for fi/faa) = ±2 separately for each of the several chosen values of 8q; (a) 
free contact; (b) pinned contact line. Note that for geometrical reasons 6i < 9q — 5 with 5 > a/Ro; 



(Fig. [T]) between the particles is almost equal to n, which, in the absence of the substrate 
(i.e., for a free droplet with a fixed center of mass), would indeed correspond to a local free 



energy minimum 



19j. This shows that the presence of the substrate, even if the particles 



are far away from the contact line, can change the interactions also qualitatively. 
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(a) free contact line, : = 2 = Q - 8 = Q - te/36 




'— 6 = 271/12 — e Q = 671/12, theory 

- 9 = 3tc/12 - = 7tc/12 

— 9 = 4tc/12 -- 9 = 8tc/12 



- 9 = 5tc/12 



9 = IOtc/12 



<t> 2 -<t>! 

(b) pinned contact line, 0=0=0 . (0„) 

v ' r 12 mm 




- O = 471/12 

- Q = 671/12, theory 

- O = 871/12 

- O = IOtc/12 



4> 2 -<K 



FIG. 8: The effective pair potentials for two particles of equal radii a, both placed at the polar 
angles 9\ = 02 corresponding to the minimum of the one-particle free energy landscape Ai 7 ^ (see 
Fig- Ej)) as a function of the relative azimuthal angle 4>2 — 4>i between the particles (see Fig. [1]) 
for various contact angles 9q. (°0 f ree contact line: 9\ = 62 = Oq — 5, with 5 = 7r/36 (Fig. 0a)); 
(b) pinned contact line: 6\ = 62 = 9 m in(Go) (Fig. 0^))- The colored lines represent the numerical 
results for a droplet of radius Ro/a = 20 for 6>o = 2-7T/T2 and #0 — 47r/12, Ro/a = 14 for 60 — Sir/ 12, 
and Ro/a = 12 for the remaining contact angles. The thick black lines correspond to the analytical 
expression V^(cj>2 — 4>i) — V (T (Q) with V c given by Eq. (fTTj) . 



19 



IV. SUMMARY AND DISCUSSION 



We have obtained analytical results for the free energy of two capillary monopoles at the 
surface of a sessile drop (Fig. [T]). These monopoles are generated by external forces acting in 
radial directions. The results, which have been obtained by using the method of images for 
the special case of a contact angle 6$ = n /2 but for arbitrary angular positions Qi = (9i, (pi) 
of both monopoles i — 1,2 (Figs. [3]and|3]), show that the preferred configuration of the 
particles depends sensitively on the boundary conditions at the substrate (see Fig. [2]) . For a 
free contact line the particles are attracted by the contact line (see Q\ = 77r/18 in Fig. |2fa) 
and Q\ = Att/9 in Fig. E](a)). In this case two locally stable configurations are possible: the 
one in which the particles touch each other and the one with the particles at opposite sides 
of the contact line. On the contrary, for a pinned contact line the particles drift towards 
a characteristic polar angle 9 min ~ 49° (see 9i = 6 min in Fig. WLp)) and there they always 
attract each other; in this case there is no other locally stable configuration. For 6q = n/2 
the effective pair potential between a probe particle 2 and a fixed reference particle 1 is 
illustrated for various reference configurations Qi in Figs. [3] and H] for a free and a pinned 
contact line, respectively. 

We have also obtained numerical results for an arbitrary contact angle 0q ^ 7r/2. For 
a single particle i the results are qualitatively the same as for 6*0 = tt/2 (compare Figs. [2] 
and [7]). There is a free energy barrier between the two minima at 9i = and 0, = 8 — 5 
for a free contact line and a single global minimum at 9i = # m m(#o) in the case of a pinned 
contact line. In the presence of two particles the corresponding effects described above for 
$o = 7r /2 are to a certain extent preserved qualitatively. However, for # > 57r/6 and with 
both particles at the free contact line, the second minimum corresponding to particles being 
at opposite sides of the contact line disappears and their mutual attraction varies purely 
monotonically (see Fig. 0(a)). 

In general the free energy barrier increases significantly with 9q (see Fig. [7(a)). For 
6*o = 5?t/6 it is one order of magnitude larger then for 6q = tt/6. The above analysis 
is valid for micrometer-sized or smaller droplets for which the effects of gravity can be 
neglected. The radial forces acting on the particles can be generated by using, e.g., optical 
tweezers. In this set-up the free energy landscapes AF& could be determined for example 
by monitoring the Brownian motion of the probe particle with the reference particle pinned 
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by an optical tweezer. Alternatively one could directly measure the tangential force on the 

(2) 

probe particle VaAi 7 ^ = const, Q 2 = £i)/Rq. The first method works successfully if 
the particles can explore all angular configurations, which is possible only if the free energy 
barrier is comparable with or smaller than the thermal energy /c#T. For particles of equal 
radii a\ = a 2 = a on a droplet with a surface tension 7 rs 0.05N/m, and subjected to radial 
forces f\ = f2 = 70, the typical free energy variations across the whole droplet are of the 
order of 7a 2 /10 (see Figs. [3]andH]), which at a room temperature correspond to variations 
of the order AF^ ~ 10 6 &bT for a = 1/im, ~ 10 4 &bT for a = lOOnm, and ~ 10 2 k bT for 
a = lOnm. These values are independent of the droplet radius Rq as long as the latter is 
smaller then the capillary length. They are reduced by two orders of magnitude if the forces 
are reduced by one order of magnitude. In such a case particles of sizes a < lOnm could 
indeed explore the whole free energy landscape. 

In view of the possibility of experimental realizations it is instructive to discuss if the 
anisotropies of the interactions induced by the presence of the substrate are relevant for 
actual systems. For Ro/a ~ 23 the polar angular separation 6*21 = 7r/36 ~ 0.087 corresponds 
to the situation that the particles almost touch each other. In the case of a pinned contact 
line and 8 = ir/2 the amplitude of the free energy variations of the probe particle as a 
function of the local azimuthal angle 02i with the reference particle positioned at 9\ ~ 8 m i n 
(see Fig. [5]) is given approximately by 10~ 3 7a 2 . For fi = f 2 = 7a this corresponds to lO^ksT 
for a = 1/im, 10 2 ksT for a = lOOnm, and &bT for a = lOnm. Therefore we expect that for 
particles of sizes a > lOnm the anisotropy of their interactions influences significantly the 
equilibrium configurations of the particles. 

Appendix A: Finite size correction to the numerically calculated free energy 

In this appendix we calculate a finite size correction to the numerically calculated free 
energy. In the reference configuration (fi = f 2 = 0) the immersed part 5V re f^ of particle 
i (being the intersection of the domain occupied by the particle with the spherical cap, 
representing the reference droplet of volume Vi = (47if (9 )Rl/3)(l+O(a/R ) 3 ) with fo(0) = 
(2 + cos6')(l — cos#) 2 /4) has, independently of Q{ and for the contact angle at the particle 
equal to 7r/2, the volume 8V re f^ = (27ra 3 /3)(l + O(a/R )). The position of this cavity in the 
liquid depends on f2j. Therefore the position xcM,ref of the center of mass of liquid depends 
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on Qi and f2 2 . In the case # 2 = Q\ it equals (here the direction of the x-axis is taken such 
that (f>2 = A0/2 and 0i = — A0/2, where A<fr = 2 — <fii) 



x 



XCM,ref(0l,4>2 ~ 4>l) = ~^-rj = ~T7 I dV X 

VI Vl JSV ref}1 +8V refa 

2SV re f t i „ . „ r/ , , , a 



3 



i2osin0icos[(0 2 -0i)/2] « - sin ^ cos[(0 2 - 00/2] . (Al) 

This gives rise to the free energy contribution SF which can be understood as the total work 
done by the force fcM{6i,<j>2 — 4>i) = — 2/i sin^i cos[(02 — 0i)/2] applied to the center of 
mass in order to counterbalance the lateral component of the forces f\ and / 2 (below we 
assume fi = / 2 ) upon displacing the center of mass by —XcM,ref(&i, 02 — <Pi)'- 

n(f>2—4>\ 

5F = d(f) fcM(<P)[-dx C M,ref(0 1 ,(f))/d(f)] 



Jia 3 



sin 2 Q x \ d(j) sin(0/2) cos(0/2) 



3 - sin 2 ^ cos 2 [(0 2 - x )/2] (A2) 



^o/o(^o) 

The results presented in Fig. [S^o) have been obtained by subtracting this correction 5F from 
the numerically calculated free energies. 
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